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Abstract 

We give a mathematically rigorous analysis which confirms the surprising 
results in a recent paper [2] of Benilov, O'Brien and Sazonov about the spec- 
trum of a highly singular non-self-adjoint operator that arises in a problem 
in fluid mechanics. 
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1 Introduction 

In a recent paper [2j Benilov, O'Brien and Sazonov have shown that the equation 

df d f . ,^,df\ Of 

approximates the evolution of a liquid film inside a rotating horizontal cylinder. 
The variable 6 is taken to lie in [— vr, vr] and one assumes that the solutions / are 
sufficiently smooth and satisfy periodic boundary conditions. 

The operator H is highly non-self-adjoint (NSA) and it is not amenable to standard 
elliptic techniques because the second order coefficient is indefinite. For 6 G (0, vr) 
the second order term has a diffusive effect on the evolution but for 6 G (— vr, 0) its 
effect is anti-diffusive. Many of the calculations in [2] are based on an asymptotic 
or WKB analysis for small e > 0, but this has dangers because infinite order ap- 
proximate eigenvalues of NSA operators need not be close to true eigenvalues. Nor 
need eigenvalues computed by truncations of a highly non-self-adjoint operator to 
large finite-dimensional subspaces by standard methods be close to the eigenvalues 
of the original operator; see [31 HJ E] for examples and discussions of their relation- 
ship to pseudospectra. Our goal in this paper is to rederive some of the results in 
[2] for a fixed positive value of e by a rigorous and non-asymptotic technique. We 
also provide strong numerical evidence that the eigenvectors do not form a basis. 
In our numerical calculations we take e = 0.1 . 



1 



Before proceeding we mention that closely related operators were discussed in [6l 
pp. 124-125, 406-408] and pLj- 



2 A Reformulation of the Problem 



We focus attention on the spectral properties of the operator 



{Hfm:=s 



d_ 




defined on all periodic functions / G 7r,7r). We normally assume that 
< e < 2 for reasons explained in Corollary [21 According to the WKB analysis of 
[2j the eigenvalue equation —iHf = Xf has a sequence of real eigenvalues which 
converge to the integers as e — > 0. This suggests that the evolution equation is 
neutrally stable, but Benilov et al. show that it exhibits explosive disturbances. 
This is closely related to the pseudospectra of the operator. Our goal is to prove 
that there are indeed real eigenvalues A without depending on WKB analysis, and 
to provide a simple and rigorous method for computing them. 

By expanding / G tt, vr) in the form 



one may rewrite the eigenvalue problem in the form Av = Xv, where A = —iH is 
given by 

{Av)n = -n{n - l)vn-i - -n{n + l)vn+i + nvn- 
The (unbounded) tridiagonal matrix A is of the form 



where A_ acts in /^(Z_), the central acts in C and Aj^ acts in /^(Z+). The 
coefficient map n — > —n induces a unitary equivalence between A_ and — so we 
only need study the spectrum of A^. Since A\ = DA^D~^ where D^,* = ^r,s{—'^Yi 
Aj^ and A\_ have the same spectrum. We assume that ^4+ has its natural maximal 
domain, and observe that it is a closed operator. We will see that its eigenvectors 
decrease more rapidly as ^ +oo the smaller e > is. We prove that the spectrum 
is discrete, i.e. that it consists only of isolated eigenvalues of finite multiplicity, in 
Section [31 

Benilov et al. correctly state in [2] that one obtains very poor numerical results 
if one simply truncates A to produce a finite matrix whose eigenvalues are then 
computed. We study the matrix A in a completely different manner. 
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V A_ 
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The eigenvalue equation for may be written in the form 

Tl — A 

n{n — l)f„_i — n{n + l)f„+i + 2 Vn = 0. (1) 

The reahty of the coefficients of ([T]) imphes that if A is an eigenvalue then so is A. 
It does not, however, imply that the eigenvalues are real. 

We confine attention to the solutions of ([T]) with support in Z"*", and regard the 
n = 1 case, namely ev2 = (1 — A)fi, as an initial condition. Since it is a second 
order recurrence equation, the solution space of ([T]) is two-dimensional. We will 
see that one solution, often called the subordinate solution, lies in P{Z^), but 
no others do so if < £ < 2. We say that A > is an eigenvalue of if the 
subordinate solution of the recurrence equation satisfies the initial condition. 

If one assumes that ([T]) has a solution of the form Vn = n°'{l + b/n + 0{l/n'^)), then 
one finds that a = —1 + 1/e and b = X/e. This motivates our next two lemmas. In 
the following calculations we introduce constants A^}*^, and will use the fact that 
they can always be increased without affecting the results. 

Lemma 1 If X > 0, there exists N = N^l such that if f„ is a solution of ^^ 
satisfying < v^^i < {N — i)"- for i = 1, 2 where a = —1 + 1/e, then < f „ < n° 
for all n > N . 

Proof Suppose that n > A + 3 and < < {n — for 2 = 1,2. Then 

„ , - .„ 1-A 

< n "u^ 



< 

= l-—^+0{n-^) 
< 1 

for all n > = A^{^]. It follows inductively that < f „ < ri" for all n> N. 

Corollary 2 // A > and e > 2 then every solution of (QP lies in Z^(Z_|_). In 
particular every such X is an eigenvalue of A^ . 

Proof Let = A^}^^. Let u be the solution of ([I]) such that un^2 = and 
un-1 = {N — 1)", and let v be the solution such that vn-2 = {N — 2)" and 
vn~i = 0. Since a < —1/2, Lemma [1] implies that both lie in P{Z^). The space of 
all solutions is two-dimensional, so every solution lies in /^(Z_(.), and this applies 
in particular to the solution that satisfies the initial condition. 

It is highly probably that one could avoid the above conclusion by imposing bound- 
ary conditions at -|-oo if e > 2, i.e. by reducing the domain of A^. We do not 
pursue this possibility. 
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(2) 

Lemma 3 For every A G R there exists N = N\ j such that if f „ is a solution of 
(Hp satisfying 

Vn>i^l + ^) (2) 

for n = N — 1 and n = N — 2, where a = —1 + 1/e and k = 1 + X/e, then the same 
inequality holds for all n > N. 

Proof Suppose that n > X + 3 and that 

Vn-i > (1 + -^) {n-iY 

for i = 1, 2. Then 

-a(n-2 n-l-X \ 

n Vn = n Vn-2 + 2 -. 7TVn~l 

\ n en{n — I) J 

> + A_^y"' + 2"7^-Vi + ^) fi-i 

y n — 2 J \ nj en[n — 1)^ n — 1 J \ n 

= l + - + 4 + 
k 

> 1 + - 

n 

provided n > N = N^l. It follows inductively that ([21) holds for all n > N. 
Theorem 4 If < e < 2 and X is a real eigenvalue of then A > 1. 

Proof Suppose that A^v = Xv where A < 1 and fi = 1. The initial condition 
ev2 = (1 — A)fi implies that f2 > 0, and it then follows from the signs of the 
coefficients in ([T]) that f „ > for all n > 3. Lemma [2] implies that there exists a 
constant c > such that 

for all n > N^^. The lower bound a > —1/2 implies that v ^ /^(Z+), and hence 
that A is not an eigenvalue of A^. 

Hypothesis From this point onwards we assume that < £ < 2 and A > 0. 

Theorem 5 For every 6 > there exists N = A^A.e,<5 (in-d a solution v of (Qp such 
that 

< Vn < (1 + 
for all n > N, where a = —1 + 1/e. 
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Proof We put N = Nx,e,6 = max{A^i^], A^'jJ, 2 + k/6} where k = 1 + \/e and let 
V be the solution of such that VN-i = {1 + S){N — iY for z = 1, 2. Lemma [H 
implies that < f„ < (1 + for all n > N. Since 

for n = N — 1 and n = N — 2, we deduce by Lemma [3] that fn > (1 + k/n)n^ > rf" 
for all n > N. This completes the proof. 

We will show that, up to a multiplicative constant, there is exactly one 'subor- 
dinate' solution t> of ([T]) such that lim„_,+oo f « = 0. We identify this solution by 
solving the recurrence relation backwards from n = M and then letting M +oo. 

Lemma 6 There exists N = N^l such that if M > N and Vn = (— is a 
solution of (QP satisfying < wm+i < {M + i)^'^ for i = 1, 2 where c = 1 + 1/e, 
then < Wn < n~'^ for all n satisfying N < n < M . 

Proof The sequence w„ satisfies the recurrence relation 

n + 2 2(n + l-A) 

Wn = Wn+2 H 7 — — m„+i . (3) 

n en[n + 1) 

This has positive coefficients for > A so the solution is positive if A < n < M. 
Suppose inductively that < Wn+2 ^ {n + 2)~'^ and < Wn+i < {n + l)~^ for such 
an n. Then 

/ 2\i-'= 2(n + 1 - A) / 1^ 
V nj en[n + Ij V n, 

= l-^ + 0(n-3) 
< 1 

for all large enough n. By induction there exists N = Nx^^ such that < w„ < 
provided N < n < M. 

Lemma 7 There exists N = N^l such that if M > N and f„ = (— l)"'tf„ is a 
solution of (QP such that 

w^>(^l- n-^ (4) 

for n = M + 1 and n = M + 2, where c = 1 + l/e and h = 1 + X/e, then ^ holds 
for all n satisfying N < n < M . 

Proof Suppose that maxj/i, X} < n < M and (jl]) holds when n is replaced by ra + l 
or n + 2. Then 

/ h \ / 2\i-^ n+1- X f h \ / 1 

n'^Wn > 1 1 + - + 2^ r 1 1 + - 

V n + 2/ V nJ enin + 1) \ n + V n 



= l-- + -,+Oin-') 
h 

> 1 - - 

n 

provided n is large enough. An induction now implies that there exists = A^}'*] 
such that (jlj) holds for all n such that N < n < M. 

Theorem 8 There exists N = N^^^l and a unique solution u„ = (— of (QP 
such that 



1 n " < Wr, < n " 



for all n > N , where c = 1 + 1/e and h = 1 + \/e. Hence 



lim Wnu'^ = 1. (5) 



n— >+oo 



Proof Let M > N = A^^J = and let w^^) denote the solution of 

such that w!^^ = for n = M + 1 and n = M + 2. Lemmas El and [7] imply 
that 

(l - n-- < wf^) < n-' 

for all n such that N < n < M. By choosing a sequence Mr —>■ +oo such that 
ly^^''^ and w^^'']^ converge as r ^ +oo we see using (|3]) that w^^^''^ converge for all 
n > 1. Denoting the limit by w^°°^ we deduce that 

1 - - 1 < < n-' 
nj - n - 

for all n > N. Putting = {—l)'^w^\ the uniqueness of the solution f*^°°) 
subject to the normalization condition ([5]) follows from the fact that the solution 
space of (lH) is two-dimensional and it contains a divergent sequence by Theorem [5l 

Numerical examples suggest that the following lemma is not the best possible 
and that w takes its maximum value very close to n = A. Figure 1 plots the 
eigenfunction v of the operator for the eigenvalue A ~ 14.94784 with e = 0.1. 
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Figure 1. Eigenvector for e = 0.1 and A ~ 14.94784 

Lemma 9 // A > then the unique subordinate solution w of ([^j satisfies < 
Wn+i < Wn for all n > 2X. 

Proof Let w^^^^ denote the solution of constructed in the proof of Theorem [HI 
Then 



(M + 1)-" f(l + 2/M)^-^(l + 1/M)^+ 



eM 



A/(M + 1)) 



= (M+l)-^(l + c/M + 0(M-2; 
> (M+1)-^ 
provided M is large enough. Therefore 



w 



(M) > (M) ^ (M) 



We prove inductively that > w^J^^l for all n such that 2X < n < M. If this 

holds with n replaced by n + 1 or by n + 2 then 



w. 



(M) 



W. 



(M) 
n+1 



n + 2 



-w 



n 

n + 3 

1 

n + 1 



(M) 
n+2 



-m 



(M) 
n+1 



-w, 



(M) 
n+3 



2(n + 1 - A) 
en{n + 1 

2(n + 2- A) 
e(n + l)(n + 2 



-w. 



(M) 
n+2 



n + 2 



(M) 



n + 3 



"+2 ^ + l^"+3 



n 

2fn+ 1 



> 



en{n + 1) 
n+2 n+3 
n n + 1 



^)„,.{M) 
"^n+1 



2(n + 2 - A) (M) 
£(n + l)(n + 2)'^"+' 



(M) 
n+3 



en{n + l) e{n + l){n + 2) 

> 

provided n > 2X. This completes the induction. 

Finally we take the same sequence M(r) as in the proof of Theorem [8] to obtain 
< wi°?] < for all n > 2A. 



3 Compactness of the Resolvent 

In this section we prove that ^ Spec(A_|_) and that A'^^ is a Hilbert-Schmidt 
operator, and hence compact. This implies that the spectrum of is discrete 
and coincides with its set of eigenvalues. We cannot, however, prove that the 
spectrum is real. We define the Hilbert-Schmidt operator R on /^(Z+) by 

oo 

{Rf)m = Pm,nfn (6) 
n=l 

where p G /^(Z+ x Z+) is given explicitly. We then show directly that R is the 
inverse of Aj^. 



Let d) be the solution of 



^-2^ 2 ^ 

\~2 H 0n-l 



n en 

that satisfies the initial conditions 0i = 1 and 02 = One sees immediately 
that 0„ > for all n > 1. Theorem O implies that there exists a constant Ci > 
such that 

Ci^n°- < (j)n < cin°- 

for all n>l. 

Let ipn = (— be the unique subordinate solution of 

n-2 , 2 , 

Wn = Vn-2 H Vn-1 

n en 

such that w satisfies the asymptotic condition hm„^+oo n^Wn = 1- Since 

n + 2 2 

Wn = Wn+2 H 

n en 

we see that w„ > for all n > 1, and indeed that there exists a constant C2 > 
such that 

C2^n~'^ < Wn < C2n~''^ 

for all n > 1. 
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We finally put 

a„ = \n{n - + |n(n + 1)0„^„^, + n0„^„ 

and observe that cr„ > for all n > 1. The upper and lower bounds on </> and w 
imply that there exists a constant C3 > such that 

for all n > 1. 

Theorem 10 If0<e<2 and 

^ J (-l)"0^^/'„/cr„ ifm<n, 
P"''^ \ (-l)>„0ja„ z/m>n. 

i/ien p G /^(Z_|_ x Z+). T/ie Hilbert- Schmidt operator R defined by ^ satisfies 
A+Rf = f for all f E 1^^+)- Indeed ^ Spec(A+) and R = A-\ 

Proof The above bounds on 0, ip, a imply that 



c^nf-n ^ if m < n, 
c^mT^n'^''^ if m > n. 



It follows that 



and then that 



3 



m=l 



CXO. 



m,n=l 



We conclude that i? is a compact operator. If {en}'^=i is the standard basis in 
/^(Z+) then a direct calculation shows that A^Rcn = e„ for all n. By using the 
fact that y4+ is closed one deduces that Ran(i?) C Dom{A+) and that A+Rf = f 
for all / G /^(Z+). We conclude from this that Ran(y4+) = /^(Z+). The bound 
< e < 2 implies that Ker(y4_|_) = {0} by Theorem HI so we finally see that 
^ Spec(A+) and that R = A'^. 



4 A-Dependence 

In this section we prove that the unique normalized subordinate solution v\^n = 
(— l)"wA,n of ([T]) provided by Theorem [H] depends continuously on A. 

We first observe that for any A > 1 the various constants N^l are uniformly 
bounded with respect to A provided < A < A. We (incorrectly) use the notation 
iv|*g to refer to the relevant upper bounds. 
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Lemma 11 //0<A<;U<A then 
for all n > A. 

Proof The positivity of w\^n for n > A foHows from the positivity of the coefficients 
of ([3]) for n > A and the positivity of wx^n for all n > = Nj^l. We only need 
only prove the central inequality above since the other two are special cases of it. 

Theorem [8] implies that if 5 > then 

Wfi,n < (1 + S)wx,n (7) 

for all n > iV = Nj^l g. This inequality persists for all n G [A, A^] by the mono- 
tonicity of the coefficients of ([H]). Since ([7]) holds for all 5 > and all n > A, the 
required inequality follows by letting 6^0. 

Lemma 12 //0<A</i<A and |yU — A| < 6 then 

< Wx,n < PA,e,n,5W^^n (8) 

for all n > 2A, where 

PA,e,n,5 = (1 + ^) GXp | 2^"^ E ^" f ' 

I. r=n ) 

Proof Since 1+5 < pA,s,n,5, TheoremOimplies that ([H]) holds for all n > A^ = N^l 5. 
We prove inductively that the same inequality persists for n G [2A, A^]. If holds 
with n replaced by n + 1 and by n + 2, then, using Lemma [HI we obtain 

n + 2 2{n + l-X) 

U)x,n = Wx,n+2 H 7 ——Wx,n+1 

n en[n + Ij 

n + 2 2(n + l-A) 

< PA,e,n+2,5Wn n+2 H 7 , n PA,e,n+l,SW n n+l 

n en[n + 1) 

(n + 2 2(n+l-A) \ 
\ n en[n + 1) j 

(n + 2 2(n+l-/i) 25 \ 

< PA,e,n+l,S W^,n+2 H 7 r~r;~^M 

,n+l H 

\ n en[n + 1) ew^ J 

< PA,e,n+l,5 ( W^,n + — W^,n+1 



/ 2(5^ 

< PA,e,n+l,5 ( 1 + — 



en^ 



< PA,e,n,SWf,^n- 

This completes the induction. 



10 



Theorem 13 The subordinate solution v\ depends continuously on A for < A < 
oo. Hence the function 

/(A) := evx,2 - (1 - A)wa,i 

is continuous on [0, oo). 

Proof It is sufficient to prove that / is continuous on [0, A] for every positive 
integer A. It follows directly from the estimates in Lemmas [TT] and [12] that the 
map A G [0,A] — (wa,2A, 'W^a,2A+i) is continuous. Composing this with the linear 
(and therefore continuous) map (u'a,2A5 'U^a,2A+i) ^v\^2 ~ (1 ~ '^)'^a,i yields the 
second statement of the theorem. 

5 Numerical Calculations 

Let v\ denote the solution of ^ such that 

lim vx,n{-irn'' = 1. 

Then A > is an eigenvalue if and only if 

/(A) := evx,2 - (1 - A)^'A,l 

vanishes. Since this function is continuous, one can compute the roots of /(A) = 
by evaluating /(A) numerically for a range of values of A. We determined the 
subordinate solution by solving ([S]), starting from M = 4000 (and also M = 8000 
to check consistency) with WM+i = {M + i)~'^ for i = 1, 2. Figure 2 plots /(A) for 
e = 0.1 and < A < 4. The eigenvalues listed in Table 1 were obtained by solving 
/(A) = numerically, and are quite close to those obtained in [2j. 



n 


A„ 


ll^nll 


1 


1.00968 


1.0189 


2 


2.07334 


1.1848 


3 


3.22978 


1.8868 


4 


4.50134 


4.3409 


5 


5.89993 


13.341 


6 


7.43194 


50.638 


7 


9.10097 


226.20 


8 


10.9092 


1152.9 


9 


12.8578 


6561.3 


10 


14.9478 


41018 


15 


27.5331 




20 


43.74 





Table 1. Eigenvalues of for e = 0.1 
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The computation is very stable and one can confidently evaluate the first ten 
eigenvalues to much higher accuracy. The list of eigenvalues found is compatible 
with the asymptotic formula A„ ~ arC where a ~ 0.53 and 7 ~ 1.44. 

However, for e = 1, the Fourier coefficients decrease much more slowly, and the 
eigenvalue calculation is correspondingly more onerous. We computed the first 
five eigenvalues for e = 1, determining the subordinate solution as before with M 
between 1000 and 32000. The apparent numbers of eigenvalues increased from 7 to 
11 as M increased in this range. For M = 4000 it appeared that the computation 
of the first five eigenvalues presented in Table 2 was reliable. 



1 1.4485 

2 4.3159 

3 8.6219 

4 14.3638 

5 21.5414 

Table 2. Eigenvalues of ^4+ for e = 1 

X 10"' 

1 r 




I I U I I I I L 

0.5 1 1.5 2 2.5 3 3.5 

Figure 2. /(A) for e = 0.1 and < A < 4 



We conclude with some comments about the conjecture in [2] that the eigenvectors 
form a basis. It seems quite plausible that they form a complete set in the sense 
that their linear span is dense. However, if they form a basis then the spectral 
projections 

Pnf = 



12 



of H must be uniformly bounded in norm, where 0„ are the eigenf unctions of H 
and 0* the corresponding eigenfunctions of H*; see [U Lemma 3.3.3]. However, it 
appears from |2j Figure 4] that the eigenfunctions 0„ concentrate more and more 
strongly around 9 = ti as n increases; the eigenfunctions 0* should concentrate 
around 6* = as n ^ oo for similar reasons. If this is indeed the case then the 
norms of the spectral projections 

IIP II ll0n|IIICll 
" l(0n,0;)| 

must diverge as n — > oo and the eigenfunctions do not form a basis. The norms of 
the first 10 spectral projections are presented in Table 1 and confirm the conjecture 
that they diverge as n increases. See |E] for another highly non-self-adjoint operator 
arising in physics for which an apparently well-behaved sequence of eigenfunctions 
do not form a basis. 
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